Sensitivity to eV-scale Neutrinos of Experiments at a 
Very Low Energy Neutrino Factory 

J.H. Cobb and CD. TunnelQ 

Subdepartment of Particle Physics, University of Oxford, 
Denys Wilkinson Building, Keble Road, Oxford, 0X1 3RH, UK 



o 

(N 

> 
O 

:z; 

00 



I 

Oh' 

CD 



> 
O 

in 



X 



A.D. Bross 

Fermilab, P.O. Box 500, Batavia, IL 60510-0500, USA 
(Dated: November 29, 2011) 

The results of LSND have yet to be confirmed at the Scr-level. An experiment is proposed utilizing 
a 3 GeV muon storage ring that would allow for both disappearance and appearance channels to 
be explored at short-baselines. The appearance channel could provide well over Scr confirmation 
or rejection of the LSND result. Other physics could also be performed at such a facility such as 
the measurement of electron-neutrino cross sections. The sensitivity of experiments at a Very Low 
Energy Neutrino Factory (VLENF) to neutrinos at the eV-scale is presented. 

PACS numbers: 14.60.Pq, 14.60.St 



Despite the experimental success of the past decade in 
establishing that neutrinos are massive and mix, whether 
they be produced in the Sun reactors [1], accelera- 
tors or the atmosphere there exist nevertheless 
some anomalies. All oscillation experiments to date can 
be explained with three neutrinos except LSND — and 
various experiments could be made to agree with either. 
The LEP collider experiments Q showed that additional 
light neutrinos cannot couple to the Z . The MiniBooNe 
experiment has been unable to deny or confirm the LSND 
result since their data that uses the same anti-neutrino 
beam polarity as LSND agrees with the background-only 
hypothesis at 0.5% and agrees with short-baseline oscil- 
lations at 8.7% 0. An experiment with a sensitivity 
greater than 5f7 is needed in order to refute or confirm 
evidence of neutrinos at the eV-scale. 

The LSND data suggests a new mass splitting Q that 
is potentially observed elsewhere. Recently, a recalcu- 
lation of reactor fluxes resulted in the reactor neutrino 
anomaly 0, Q which, along with the Gallium anomaly 
[loj . are possiljly caused by sterile neutrinos. Global fits 
reconcile these anomalies into a framework that intro- 
duces a sterile neutrino at the eV-scale, e. 



ditional indications come from cosmology where WMAP 
favors more than three neutrinos (isf . 

The Very Low Energy Neutrino Factory (VLENF) 
would utilize a muon storage ring to study eV-scale oscil- 
lation physics and measure cross sections (including Ve)- 
Pions are collected from a target and injected into a stor- 
age ring where they decay to muons. The storage ring 
is optimized for 2 GeV muons where the energy is opti- 
mized for the needs of both oscillation and cross section 
physics. The muons decay according to /i"*" — e^v^v^- 
Straight sections in the storage ring result in neutrinos 
directed at a near- and far-detector. 

The storage ring could be a fixed field alternating gra- 
dient (FFAG) lattice for large momentum acceptance, 



which is important given the large momentum spectrum 
after the target. The pricing and engineering of FFAGs 
are well-understood because of experience building them 
14| . The design would require only normal conducting 
magnets which simplifies construction, commissioning, 
and operations. Design work for the injection into the 
storage ring and the particle collection downstream of 
the target are underway. 

The near detector will be placed at 20-50 meters from 
the end of the straight and will measure neutrino-nucleon 
cross sections of interest to future long-baseline exper- 
iments including the first precision measurement of 
cross sections. The far detector at 800 m would measure 
disappearance and wrong-sign muon appearance chan- 
nels. The detector would need to be magnetized for the 
wrong-sign muon appearance channel. Numerous possi- 
bilities exist for detector technologies that include liquid 
argon, MINOS inspired, and totally active scintillating 
detectors. For the purposes of this study, a detector in- 
spired by MINOS but with thinner plates is assumed. 
The experiment will take advantage of the "golden chan- 

y,, where the result- 
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ing final state is a wrong-sign muon. 

The probability v^. — )> depends on the mixing ma- 
trix, U . Let Rij be a rotation between the ith and jth 
mass eigenstates without CP violation. For N neutri- 
nos, Rij has dimension N x N. By convention, the three 
neutrino mixing matrix is C/pMNS — R2d,Rij,Ri2- In the 
(3+1) model of neutrino oscillations, extra rotations can 
be introduced such that the mixing matrix is J7(3-)-i) = 
i?34i?24-Ri4C^PMNS- Given that Am|j >> ArTig^, J/pmns 
can be approximated by the identity matrix (ie. the 
"short-baseline approximation"). It then follows that 
|J7e4p — sin(0i4) and |?7^4p = sin(6'24) cos(0i4). 

The oscillation probabilities for appearance and disap- 
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pearance, respectively, are: 



= 1- [4|[/a4p(l-|C/a4p)] sin' 
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Disappearance measurements will constrain |C/;i4p and 
|J7e4p while the appearance channel measures their prod- 
uct |J7e4p|^^/i4p- These parameters are over-constrained. 
Only the appearance channel will be assumed since 
the physics sensitivity of disappearance measurements is 
well-understood for short-baseline experiments. Without 
oscillations there would be ~ 10^ charge current (CC) 
events in the far detector. The disappearance oscillation 
probability, if LSND is correct, would be on the order of 
a few percent for this experiment; the statistical errors 
and systematic errors are expected to be about half a 
percent from the presence of a near detector and beam 
instrumentation. 

The oscillation sensitivities have been computed us- 
ing the GLoBES software (version 3.1.10) [ll[l3|. Since 
GLoBES, by default, only allows for a 3 x 3 mixing ma- 
trix, the SNU (version 1.1) add-on [13, [11] is used to 
extend computations in GLoBES to 4 x 4 mixing matri- 
ces. 

Unlike previous analyses, the far detector approxima- 
tion of the source and detector being treated as point- 
sources cannot be used since the size of the detector 
and accelerator straight are comparable to the baseline. 
This study computes the neutrino flux by integrating the 
phase space of the stored muons using Monte Carlo (MC) 
integration. The beam occupies a 6D phase space {x, y, z, 
Px,Py, Pz) and the detector has a 6 m x 6 m cross section. 
A random point is chosen within the beam phase space 
and within the detector volume. The transverse phase 
space is represented by the Twiss parameters a = and 
/3 = 25 m where the la Gaussian emittance is assumed 
to be 15 mm. It follows that the spread in, for exam- 
ple, X is ax — and the angular divergence in x is 
fa;' = \ftf^. The longitudinal phase space (z and Pz) 
is described by assuming a uniform distribution in z and 

= (2 ± 20%) GeV. The code for the analysis herein is 
available online [3| under the GPL license [20| . 

This analysis assumes 2 x 10^^ decays of pf^ in each 
straight, normalized to an exposure similar to MiniBooNe 
of 10^^ protons on target (POT). The far detector at 
800 m consists of a kilotonne of fiducial target mass. 
The background rejection for the wrong-sign muons is 
assumed to be 10~^ for the — v^^ neutral current 
(NC) events and 10~^ for the charge misidentification of 
v^l — ^ v^;^ CC events. The backgrounds from v^. — > Vf, CC 
and NC are negligible. The v,, — >■ z/^ CC events have a 
90% detection efHciency. These efficiencies were extrap- 
olated from previous studies 2]| and reconfirmation of 
these calculations for energies below 2 GeV is still re- 
quired. 
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FIG. 1. The oscillation probability for the "golden channel" 
— ^ i^/j using the (3-1-1) oscillation parameters in TABLED 
A baseline of 800 meters is assumed. This numerical treat- 
ment of the binned oscillation probability agrees with Eq. [1] 
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FIG. 2. The spectrum of signal and background events at the 
far-detector. The number of wrong-sign muon appearance 
signal events is 27 which arises through the "golden channel" 
oscillation probability v^j,. Roughly 2 background events 

are in the signal box and correspond to both Vy, CC with 
charge misidentification and NC events where a short-lived 
pion has faked a wrong-sign muon. An exposure of 10^^ POT 
is assumed. 



To initially gauge the sensitivity of this experimental 
setup, some oscillation parameters must be assumed that 
include the LSND effect. The best-fit data for (3-hl) 
sterile neutrinos is used [2^ where the MB v and LSND v 
data was fit (TABLE [l| . Lower energy neutrinos contain 
the most information about the mixing matrix (FIG.[T]). 

The best-fit values allow for an estimation of the event 
rates. The true number of events without efficiencies is 



TABLE I. Best-fit oscillation parameters for the (3-1-1) sterile 
neutrino scenario for MB v and LSND v data [22|. 



Parameter 


Value 


Am^i [eV"] 


0.89 


|t/e4p 


0.025 
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FIG. 3. Optimization of the mean stored muon energy and 
far detector baseline. Above muon energies of 1 GeV, appear- 
ance physics is not sensitive to the muon energy. For a fixed 
neutrino energy, the boosted muon frame will result in more 
events the more it is boosted. The optimal baseline is around 
800 meters. 

shown in TABLE |IT1 Applying the assumed efficiencies 
of 10-5 and 10"'' for ^ CC and Dx -> i>x NC, 
respectively, reveals the event spectrum (FIG. [2]). After 
cuts, there are 27 signal events and 2 background events. 

In practice the muon energy is fixed by the needs of 
the cross section physics. However, given a fixed fiducial 
mass, the baseline and muon energy can be optimized 
for short-baseline physics. The metric for comparison is 
the between simulations for the (3-1-1) best- fit men- 
tioned earlier and the background-only hypothesis. The 
short-baseline oscillation physics reach is comparable at 
energies above 2 GeV and the optimal baseline is around 
800 meters (FIG.E]). 

Having chosen a baseline and energy, it is possible to 
investigate the physics reach of this facility if we ap- 
propriately define the desired value. The design re- 
quirement is to measure the LSND effect at "5cr" but 
that can be ambiguously defined. Herein "ncr" is de- 
fined according to the Gaussian probability of the re- 
spective deviation. Let CDFx be the cumulative distri- 
bution function of some distribution X. If G is a Gaus- 
sian with /i = and cr = 1, then the p- value of a 5ct 
effect is 1 - CDFg(5) ~ 3 x IQ-^. Similarly, if xl is 
the chi-squared distribution with two degrees of freedom. 



TABLE II. A table of the raw event rates. The first row corre- 
sponds to the "golden channel" appearance signal. The other 
rows are potential backgrounds to the signal. The background 
that drives the analysis is — > S'fi. 



Channel 


Interaction 


Pre-cuts 


Ve V^l 


CC 


30 


i>x — ^ vx 


NC 


16850 


V^i 


CC 


42545 


Ve Ve 


CC 


78974 




lO'' 10"^ 10^^ 10^1 

sin' (26,^) 

FIG. 4. Sensitivity to sterile parameters. The orange band 
corresponds to the LSND v and MiniBooNe v sterile neutrino 
fit performed by Giunti and Laveder [lH where they assumed 
only one extra massive neutrino. The contours are the sensi- 
tivities of an experiment at a 3 GeV FFAG storage ring using 
parameters defined in the text. The entire 99% confidence 
interval for the LSND v and MiniBooNe v results would be 
confirmed or excluded at 1g. 



then the required value of the x^ to be an no effect is 
CDF ,i(CDFG(n)). 

The x^ is computed using the pull-method, which al- 
lows for systematics to be included. The signal and back- 
ground normalization errors are assigned to be 2% and 
20%, respectively, and these systematic uncertainties are 
marginalized over. Spectral information is also used. The 
sensitivity to new physics is shown in FIG. |H The 99% 
confidence interval of a fit in a (3-1-1) scheme to the Mini- 
BooNe V and LSND v data is shown and excluded at la. 
This meets and exceeds the requested criterion for future 
eV-scale experiments. 

Some effects are ignored that could reduce the signifi- 
cance of the result. For example, not included are back- 
grounds associated with cosmic muons. This could be 
problematic since muons fill the storage ring which re- 
sults in a large duty factor. One solution, if it is an issue, 
would be to put RF cavities in the storage ring to bunch 
the beam. This merits further study. 

It has been demonstrated that a VLENF has sensitiv- 
ity to eV-scale oscillations. Further optimizations of the 
target horn, collection, injection, and the dynamic aper- 
ture of the storage ring will increase the neutrino flux. 
But previously unthought-of backgrounds may arise and 
detector charge identification performance is difficult at 
1 GeV and below. A design study must be performed to 
ensure enough contingency in the sensitivity to guaran- 
tee a firm confirmation or refutation of the LSND effect. 
Further study is required. 

Thanks for guidance and useful discussions with 
Joachim Kopp, David Neuffer, Kenneth Long, and Alain 
Blondel. Also thanks to Nick Ryder for a careful reading. 



4 



* [Correspond ing auth or: tunnell@fiial.gov| 
[1] B. Aharmim et al. 



Phys. Rev. Lett. 101, 111301 (2008) 
arXiv:0806.0989 [nucl-ex' 



(SNO), 



[2] M. H. Ahn et al. (K2K), Phys. Rev. D74, 072003 (2006) | 

,arXiv:hep-ex/0606032 . 
[3] R Adamson et al. (MINOS), 



[4] 

[5] 
[6] 

[7] 

[8] 



Phys. Rev. Lett . 101, 131802 (2008) 
arXiv:0806.2237 [hep-e^ 

R. Wendell et al. (Kamiokande), 

Phys. Rev. D81, 092004 (20l0)] 
arXiv;1002.3471 [hep-ex], 

Phys. Rept. 427, 257 (2006)] |arXiv:hep-ex/0509008 



A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration) , 
|Phys. Rev. Lett . 105, 181801 (201 0) 
C. Athanassopoulos ei al. (TTSND Collaboration), 

[Phys. Rev. Lett. ^ 1774 (1998) 

T. A. Mueller al, Phys. Rev. C83, 054615 (2011) 



[9] R 



|arXiv:110O663 [hep-ex 



Huber, 



Phys. Rev. C84, 024617 (2011) 



arX iv:1106.0687 [hep^ 

[10] M. A. Acero, C. Giunti, and 

M. Laveder, | Phys. Rev. D 78, 073009 (2008)| 
[arXiv :071 1.4222 [hep-phjl 
[11] C. Giunti and M . Laveder, (2011), 
[arXiv: 1107. 1452 [hepiph]] 



[12; 

[13; 

[14 

[is; 

[16; 

[17 

[is; 

[19 



[20 
[21 

[22] 



J. Kopp, M. Maltoni, and 
Phys. Rev. Lett. 107, 091801 (2011) | 
arXiv: 1103.4570 [hep-phj] 



E. Komatsu et al. 

Astrophys. J. Suppl. 192, 18 (2011)[ 



T. Schwetz, 



(WMAP), 



M. Rolinec, and 



'arXiv:1001.4538 [astro-ph.CO] 

http : //hadron . kek . j p/FFAG/FFAG05_HP/index . htm I 
(2005), see talks including the summary talk. 

P. Huber, M. Lindner, and W. Win 

ter, [C omput.Phys.C ommun. 167, 195 (2005) 

arXiv:hep-ph/0407333 [hep-ph] 
P. Huber, J. Kopp, M. Lindner, 

W. Winter, Comput.Phys.C ommun. 177, 432 (2007)| 
arXiv:hep-ph /070118 7 [hep-ph]' 

J. Kopp, I nt. J. Mod . Phys. C19, 523 (2008), erratum 
ibid. C19 (2008) 845, arX iv:physics/0610206, 
J. Kopp, M. Lindner, T\ Ota, and 

Sato, P hys. Rev. D77, 013007 (2008) 

|arXiv:0708.0152 [hep-ph] 

C. D. Tunnell, https : //code . launchpad. net/- c-tunnelll/+junk/v] 

(2011), questions should be directed to the corresponding 

author. 



http : / /www . gnu . org/licenses/gpl-3 . . html 

A. B. Laing, "Optimisation of detectors for the golden 
channel at a neutrino factory," (2010). 
C. Giunti and M. Laveder, (2011), 
hep-piill 



arXiv: 1109.4033 



